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INTRODUCTION 


The  modifications  of  surface  waves  by  a  depth-varying  current 
field  are  described  mathematically  by  the  inviscid  Orr-Sommerfeld 
equations  of  hydrodynamic  stability  (Peregrine  (3)).  The  solution  of 
these  equations  is  required  to  derive  the  wave  parameters  of  interest 
(the  amplitude,  frequency,  direction,  and  wavelength  of  the  wave). 

This  solution  can  be  difficult  to  obtain,  even  numerically. 

An  alternative  approach  to  the  analysis  of  the  modifications  of 
surface  waves  by  depth-varying  current  fields  is  developed  here  .  The 
approach  is  based  on  an  approximate  dispersion  relation  for  wave- 
current  interactions  derived  from  the  governing  Orr-Sommerfeld  equations 
(Skop  (5))  coupled  with  the  wave  kineraatic/wave  action  formulation  of 
surface  wave  propagation  (Bretherton  and  Garrett  (1),  Crapper  (2)). 

The  approach  is  particularly  useful  in  that  it  focuses  directly  on  the 
wave  parameters  of  interest  and  eliminates  the  requirement  to  solve 
the  Inviscid  Orr-Sommerfeld  equations  to  derive  these  parameters . 

The  validity  of  the  approach  is  demonstrated  by  comparisons  with 
exact  Orr-Sommerfeld  solutions  . 

THE  GOVERNING  EQUATIONS 

Let  t  deonote  time  and  let  x^ ,  X2 ,  and  z  define  a  Cartesian 
coordinate  system  with  z  *  0  as  the  undisturbed  free  surface  of  the  water 
column  and  decreasing  with  depth.  Within  some  region  of  this  space,  a 
current  field  is  imposed.  The  horizontal  components  of  the  current 
field  are  denoted  by  Uj  and  U2  and  the  vertical  component  by  W .  The 
current  field  is  allowed  to  vary  in  space  and  time. 

A  surface  wave  is  taken  to  propagate  from  some  region  of  space 
where  the  current  field  vanishes  into  the  region  of  space  where  the 
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current  field  is  imposed.  We  wish  to  ascertain  how  the  surface 
wave  is  modified  by  the  current  field. 

We  assume  only  that  the  imposed  current  field  varies  slowly  in  the 
horizontal  coordinates  and  time  relative  to  the  wavelength  and  period 
of  the  surface  wave.  Then,  as  viewed  by  the  surface  wave,  the  imposed 
current  field  appears  to  vary  only  in  the  z  coordinate  over  at  least 
several  wavelengths  of  propagation  or  periods  of  oscillation  of  the 
surface  wave.  Further,  under  the  slowly  varying  assumption,  the 
continuity  equation  for  the  imposed  current  field  appears  locally  as 
3W/3z  *  0  to  the  surface  wave.  Hence,  as  viewed  by  the  surface 
wave,  W  can  be  put  equal  to  zero  without  loss  of  generality. 

We  denote  the  velocitites  and  pressures  associated  with  the  surface 
wave  as  uf,  uj,  w* ,  and  p*  and  write  the  total  fluid  velocities  and 
pressure  as  the  sum  of  those  quantities  associated  with  the  surface 
wave  and  those  quantities  associated  with  the  imposed  current  field. 
Substituting  the  total  fluid  velocities  and  pressure  into  the  Euler 
equations,  linearizing  with  respect  to  the  surface  wave  quantities, 
and  applying  the  slowly  varying  assumption  to  the  imposed  current  field, 
we  find  the  equations  that  govern  the  local  variations  of  the  surface 
wave  as 


3u* 

.  *  .  * 
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dw  dw  dW 
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3t  3xi  3x2 


1 

p  3z 


(Id) 


where  p  specifies  the  density  of  the  water.  At  the  mean  free 

surface  z  =  0 ,  the  Local  kinematic  and  dynamic  boundary  conditions  on  the 

surface  wave  are  obtained  as 

3n*  3n*  3n*  * 

- +  U, - +  U9 - w  -  0  (2a) 

3t  1  3xx  1  3x2 


p  -  Pgn  +  Y 


32n* 


32n* 


(2b) 


Here,  n*  is  the  elevation  of  the  surface  wave,  g  is  the  acceleration  due 
to  gravity,  and  y  is  the  coefficient  of  surface  tension. 

We  seek  wavelike  solutions  to  equations  (1)  and  (2)  as 


* 

U1 

=  u^(z)  explKk^Xj  +  k2x2  ~ 

<ot)  1 

(3a) 

* 

u2 

=  u2(z)  exp[i(k^Xi  +  k2x2  ~ 

cot)) 

(3b) 

* 

w 

*  w(z)  exp[i(k^xi  +  k2x2  “ 

cot) ) 

(3c) 

P* 

=  p(z)  exp[i(k1x1  +  k2x2  “ 

iot)) 

(3d) 

* 

n 

=  n  exptKk^x^  +  k2x2  -  <*>t) 

) 

(3e) 

where  kj  and  k2  are  the  components  of  the  wave  number  vector  of  the  surface 
wave  and  where  co  is  its  radial  frequency.  In  accord  with  the  slowly 
varying  nature  of  the  imposed  current  field  relative  to  the  surface  wave,  k.]_, 
k2,  co,  and  the  amplitude  functions  of  equations  (3)  can  vary  with  xj,  x2, 
and  t  but,  again,  only  slowly  over  at  least  several  wavelengths  of  propaga¬ 
tion  or  periods  of  oscillation  of  the  surface  wave.  Substituting  equations 
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approach  to  the  analysis  of  the  modification  of  surface  waves  by  depth- 


varying  current  fields  .  The  approach  is  particularly  useful  in  that  it 


focuses  directly  on  the  wave  parameters  of  interest  and  eliminates  the 


requirement  to  solve  the  inviscid  Orr-Sommerfeld  equations  to  derive  these 


parameters  . 


AN  APPROACH  TO  THE  ANALYSIS  OF  THE  MODIFICATION  OF  SURFACE  WAVES  BY 


DEPTH-VARYING  CURRENT  FIELDS 


The  Wave  Kinematic  Equations 


The  surface  wave  phase  function  $  is  defined  by 


<t>  =  kix^  +  k2x2  ”  mt 


Recalling  that  kj,  k2,  and  u)  are  slowly  varying  funcitons  of  x1,X2»  and 


t ,  we  have 


3  <p  dp  dip 

k  i  = - ,  k2  = - ,  oj  =  -  — 

3xj  3x2 


Taking  the  mixed  second  derivatives  of  <p,  we  find  the  wave  kinematic 


equations  as 


3k ^  3k 2 
3x2  ^X1 


3k^  dw 
—  +  —  *  0 
3t  3x^ 


3k2  3m 

- +. - -  0 

3t  3x2 


r  y  % 


Equations  (8)  are  not  linearly  independent  and  serve  to  determine  only  two 
of  the  three  phase  measures  .  The  dispersion  relation  for  the  surface 
wave  is  required  to  complete  the  wave  kinematic  equations. 


An  Approximate  Dispersion  Relation  for  Wave-Current  Interactions 

Following  some  algebra,  equations  (4)  can  be  reduced  to  a  single 
equation  for  w  as 


32Ui 


32w  2 

-  (k2  + 

3z2 


Q  u2 

k  _  +  k2— -- 


3z 


3z^ 


)  w  =  0 


kfUi  +  k2U2  ~  u 


(9) 


where  we  have  denoted  the  absolute  value  of  the  wave  number  vector  by 
7  7  1/7 

k  =*  (k£  +  k2)A/  .  Similarly,  equations  (5)  can  be  reduced  to  a  single 
condition  on  w  at  the  mean  free  surface  z  =  0  as 


(k^  +  k2U2  -  ai)2 


3w 

3z 


yk2  3U1  3U2 

[k2(g  + - )  +  (kill,  +  k2U2  -  a>) (kj^ - +  k2 - )]w 

p  3z  3z 


(10) 


Equations  (9)  and  (10)  can  be  further  simplified  by  introducing  the 


effective  current  field  U  defined  by 


kl  k2 

u  =  —  ux  +  —  U2  (11) 

k  k 


Denoting  the  surface  wave  celerity  by  c  =  o»/k ,  we  then  find 
32U 

32w  8z2 

- (k2  + - )  w  »  0  (12a) 

3z2  U  -  c 


while,  at  the  free  surface 


0  3w  yk2  3U 

(U  -  c)2  —  -  ( (g  +  - )  +  (U  -  c)  — ]  w 

3z  p  3z 


(12b) 


Let  us  assume,  for  the  moment,  that  the  magnitude  of  the  effective 
current  field  U  is  smaller  than  some  characteristic  celerity  cq  of  the 
surface  wave  .  We  can  then  write 


U  =  e  co  U 

where  U  is  a  dimensionless  velocity  of  order  unity  and  where  e  is  a 
dimensionless  smallness  parameter.  Also,  we  can  seek  perturbation 
solutions  for  c  and  w  as 


(13a) 


C  =  CQ  +  +  ... 


(13b) 


w  =  wq  +  ewj  +  ... 


(13c) 


Substituting  equations  (13)  into  equations  (12)  and  equating  like  powers 
of  e,  we  find  the  boundary  value  problem  for  wq  as 


- -  kzwA  =>  0 


(14a) 


3wq  Y* 

Cq - (g  + - )wq  =  0  at  z  -  0 

3z  p 


(14b) 


and  that  for  wj  as 


- -  kw,  »  - - wn 

3z2  1  3z2  ° 


(15a) 


2  ^1  Yk 
C  3z~  ’  (*  +  T>W1 


2c0(c0U  ■  cl)  g"  'c0  j;w0  at  z  -  0 


(15b) 


7 


The  solution  to  equation  (14a)  that  vanishes  as  z  -00  is 


(16a) 


and,  from  equation  (14b),  we  find 


eg  «  cr(k)/k 


(16b) 


where 


cr(k)  =  /gk  (1  + - )l/2 

Pg 


We  recognize  o(k)  as  the  dispersion  relation  for  the  surface  wave 


in  a  region  of  space  where  the  current  field  vanishes  and  eg  as  the 


surface  wave  celerity  in  this  same  region  of  space.  Further,  we  see 


that  our  assumption  that  the  effective  current  field  U  is  smaller  than 


cq  is  satisfied  for  very  long  deepwater  waves  (k  ♦  0)  and  for  very 


short  capillary  waves  (k  -*■<*)  .  If  surface  tension  is  neglected,  the 


assumption  is  still  satisfied  for  very  long  deepwater  waves 


Substituting  for  wg  In  equation  (15a)  and  using  variation  of 


parameters,  we  obtain  the  particular  solution  for  w^  as 


w  =  -  \  1-  e2k?  d? 

1  _•  n 


(18a) 


and,  from  equation  (15b),  we  find,  after  some  algebraic  manipulation  and 


integration  by  parts, 


-  2k  cQ  /  U  e2kzd 


(18b) 


Substituting  equations  (16b)  and  (18b)  into  equation  (13b)  and  using 
equation  (13a),  the  perturbation  solution  for  the  surface  wave  celerity 
is  determined  as 


c 


o(k)  0 

- +  2k  f  U  e^kz  dz 

k  -oc 


(19) 


and  replacing  the  effective  current  field  U  by  its  definition  from  equation 
(11),  we  find  the  approximate  dispersion  relation 


a)  =  o(k)  +  2k  [k. 


0 


(20) 


Equation  (20)  was  first  derived  by  Stewart  and  Joy  (6)  and  later  modified  to 

account  for  water  of  finite  depth  by  Skop  (3).  Skop  also  noticed  that 

equation  (20)  was  not  limited  to  very  long  deepwater  waves  or  very  short 

capillary  waves.  For,  on  integrating  equation  (20)  repeatedly  by  parts, 

one  obtains  j 

ki  V-2  302s  1 

to  =  a(k)  +  (kiU^g  +  k2U2g)  -  ( — - +  — - )  +  0(-)  (21) 

2k  9z  2k  3z  k 

where  the  subscript  "S"  denotes  the  current  and  its  derivatives  at  the  free 
surface.  Equation  (21)  is,  through  terms  in  the  surface  shear,  identical 
to  the  asymptotic  dispersion  relation  obtained  by  Peregrine  and  Smith  (4)  ' 

for  short  gravity  or  gravity-capillary  waves  (k  »  1)  riding  on  a  depth-  j 

varying  current  field.  Hence,  equation  (20)  provides  a  doubly  asymptotic  j 

approximation  to  the  exact  dispersion  relation  for  surface  waves  on  a  | 

depth-varying  current  field.  (( 

Given  a  specified  current  field,  equations  (8)  and  (20)  serve  to 

determine  the  three  phase  measures  —  ki,  k2  and  to  —  for  the  surface  i 

<i 

wave  throughout  xi ,  X2  and  t.  To  complete  the  analysis  of  the  modifications  • 

i 


of  the  surface  wave  by  a  depth-varying  current  field,  an  expression  allowing 
the  calculation  of  the  surface  wave  ampH *ude  n  is  required. 


(22a) 


(23a) 


The  wave  action  obeys  the  conservation  law 


9A  a 
3t  3x^ 


(cg^A) 


a 

+ - (c„2A)  =  0 

3x2  8 


(24) 


Here,  Cg^  and  cg2  are  the  components  of  the  group  velocity  of  the  surface 
wave  and  are  found  from  the  dispersion  relation  as 


3o) 


3o) 


(25) 


Equation  (24)  was  originally  derived  by  Bretherton  and  Garrett  (1)  on 
the  basis  of  Hamiltonian  dynamics  .  A  more  physically  motivated  derivation 
of  this  equation  can  be  found  in  Crapper  (2). 

Equation  (24)  provides  the  required  expression  for  determining  the 
surface  wave  amplitude  n  throughout  x^ ,  X2  and  t . 

Remarks 

Equations  (8)  are  purely  kineraatical  in  nature  whereas  equation  (24)  is 
purely  dynamical  .  The  accuracy  of  the  approach  that  has  been  developed 


here  for  the  analysis  of  the  modifications  of  surface  waves  by  depth- 
varying  current  fields  thus  depends  entirely  on  the  accuracy  of  the 
approximate  dispersion  relation  defined  by  equation  (20).  This  disper¬ 
sion  relation  appears  in  two  contexts  in  the  approach:  first,  to  complete 
the  wave  kinematic  equations  and  second ,  to  determine  the  components  of 
the  group  velocity  that  arise  in  the  wave  action  equation  . 

As  has  already  been  noted,  equation  (20)  is  asymptotic  to  the  exact 
dispersion  relation  for  waves  on  a  depth-varying  current  field  for  long 
deepwater  waves  (k  -*■  0)  and  for  short  gravity  or  gravity-capillary  waves 
(k  »  1).  Hence,  the  overall  accuracy  of  the  approximate  dispersion 
relation  is  determined  by  its  ability  to  mimic  the  exact  dispersion  relation 
at  intermediate  wave  numbe,  i  .  Skop  (5)  has  demonstrated,  by  example,  that 
this  ability  is  excellent  for  two  depth-varying  current  fields  for  which 
Taylor  (7)  obtained  exact  dispersion  relations  in  his  study  of  hydraulic 
breakwaters  .  The  first  of  these  was  a  uniform  current  extending  from 
the  surface  to  some  depth  d,  while  the  second  was  a  uniformly  sheared 
current  extending  from  the  surface  and  vanishing  at  some  depth  d  . 

We  wish  to  examine  here  a  significantly  more  complex  depth-varying 
current  field  than  either  of  the  two  current  fields  previously  considered 
by  Skop  . 

THE  SUBMERGED  JET 

We  consider  the  depth-varying  current  field  defined  by 


•VI 
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Uo 

- [a  +  -) 

l-o  d 


Uo  2 

- [2  -  a  +  -] 

1-a  d 


0  >  z  >  -  ad 


-ad  >  z  >  -d 


(26a) 


-d  >  z  >  -(2  -  a)d 


-(2  -  a)  d  >  z 


U2  =  0 


(26b) 


As  shown  In  Figure  1,  this  current  field  represents  a  symmetric  submerged 
jet  .  The  maximum  velocity  in  the  jet  is  -Uq  and  occurs  at  z  »  -d  . 

The  width  of  the  jet  is  2d(l  -  a)  and  the  jet  velocity  goes  to  zero 
at  z  3  -ad  and  z  *  -2(2  -  a)d .  If  one  desires,  the  jet  can  be 
thought  of  as  a  model  of  an  internal  wave  trapped  at  a  thermocllne  . 


-  (2-cr)d 


Figure  1.  Schematic  of  the  current  field  characterizing  a  symmetric 
submerged  jet. 


. 


We  take  a  surface  wave  propagating  from  x  •  into  the  jet. 

Then,  In  the  surface  wave  phase  function,  we  can  put  kj  *  k  and  k2  *  0 
without  loss  of  generality.  We  also  assume,  for  this  example,  that 
surface  tension  is  unimportant  and  set  the  coefficient  of  surface 
tension  y  =  0.  Substituting  equations  (17)  and  (26)  into  equation 
(20),  we  obtain  the  approximate  dispersion  relation  as 

p 

a  =  /k - e-2a<[l  -  +  e"4(1“a)<] 

2(1  -  a)  (27) 

where  the  dimensionless  frequency  Q,  dimensionless  wave  number  ic, 
and  the  Froude  number  F  are  defined  by 

0  =  /d7g  ui,  k  =  kd,  F  =  U0//id  (28) 

The  exact  dispersion  relation  for  the  submerged  jet  is  derived  in  the 
Appendix.  We  find  it  satisfies  the  fifth  order  polynomial  in  Q 

p5o5  +  p4n4  +  p3n3  +  p2n2  +  Pln  +  p0  -  o  (29) 

where 

p5  =  4(1  -  a)3  (30a) 

P4  «  2(1  -  ct)2F  {[e-2(1“a),c  -  2)e~2<  +  2(1  -  a)x  +  e“2aK} 

( 30b) 

p3  -  (1  -  a) { 2 F2 [ ( 1  -  a)<  -  1  +  e-2(l-<*)<) 

+  2F2( (1  -  a)<e-2<1-'a><  +  e-2*1-®)*  -  He”2* 

-  F2(l  +  e-2aK)(e-2(1_a)<  -  2)  e_2(1_a)K 
+  F2 ( 2( 1  -  a)<  -  1] (1  +  e-2aK)  -  4(1  -  a)2*} 


(30c) 


p2  =  F { 2 ( 1  -  a)2K[e“2^1_ot^l<:  -  2]  e“2ic 

-  2(1  -  o)2<[2(1  -  a) k  -  1] 

-  F2(l  +  e“2a<)  ((1  -  a)<  e-2a~a)K 

+  e-2(l-a)K  _  i]e-2(l-a)i< 

+  F2(l  +  e-2a<)  [(i  -  a)ie  -  1  +  e"2(1-a)K] 

+  2(1  -  a)2  K(e“2ctK  -  1)}  ( 30d) 

Pl  -  (1  -  a)F2{2ict(l-a)tce~2(1"a)K: 

+  e-2(l-a)K  _  i]e-2< 

-  ic(e“2aK  -  l)[e“2(l“a)<  -  2]e"2(1“°)|c 

-  2k[ (1  -  o)k  -  1  +  e"2(1"a)'c]e"2o'c 

+  k[2(1  -  a)  k  -  1]  (e“2aK  -  1)}  (30e) 

p0  =  F3K(e“2a<  -  1)[(1  -  cOice-^1*"01^* 

+  e"A<1“«)K  +  (1  -  a)<  -  1]  ( 30f ) 


As  mentioned  previously,  the  roots  of  equation  (29)  as  a  function  of 
k  for  fixed  values  of  a  and  F  provide  information  on  both  the 
dispersion  relation  for  the  surface  wave  and  the  stability  of  the  imposed 
current  field  .  The  root  that  is  relevant  as  being  the  dispersion  relation 
for  the  incoming  surface  wave  is  identified  by  its  asymptotic  behavior 
from  equation  (21)  that  Q  +  as 

Figures  2,  3,  4,  and  5  show  comparisons  of  the  approximate  and  exact 
dispersion  relations  for  the  symmetric  submerged  jet.  In  Figure  2,  the 
comparisons  are  made  for  a  =  0  and  values  of  F  -  0.25,  0.50,  and  1.00 
which  represent  a  progressively  stronger  jet.  In  Figures  3,  4,  and  5,  the 
value  of  F  is  fixed  (at  0.25,  0.50,  and  1.00,  respectively)  and  a  takes 
the  values  a  ■  0.00,  0.05,  and  0.20  which  represent  a  progressively 


Figure  2.  The  dispersion  relation  for  water  waves  on  a  symmetric 
submerged  jet  as  a  function  of  the  Froude  number  F  with  the  submergence 
parameter  a  =  0.  Increasing  values  of  F  correspond  to  progressively 
stronger  jets.  Exact  solution  ( - );  approximate  solution  ( - ). 


Figure  3.  The  dispersion  relation  for  water  waves  on  a  symmetric  sub¬ 
merged  jet  as  a  function  of  the  submergence  parameter  a  with  the  Froude 
number  F  =  0.25.  Increasing  values  of  a  correspond  to  progressively 
narrower  jets.  Exact  solution  ( - );  approximate  solution  ( - ). 


Figure  4.  The  dispersion  relation  for  water  waves  on  a  symmetric  sub¬ 
merged  jet  as  a  function  of  the  submergence  parameter  a  with  the  Froude 
number  F  =  0.50.  Increasing  values  of  a  correspond  to  progressively 
narrower  jets.  Exact  solution  ( - )  ;  approximate  solution  ( - ). 


Figure  5.  The  dispersion  relation  for  water  waves  on  a  symmetric  sub¬ 
merged  jet  as  a  function  of  the  submergence  parameter  a  with  the  Froude 
number  F  =  1.00.  Increasing  values  of  a  correspond  to  progressively 
narrower  jets.  Exact  solution  ( - );  approximate  solution  ( - ). 


narrower  jet.  In  each  figure,  the  dispersion  relation  £2  ■  Sk, 
applicable  to  a  region  where  the  current  field  vanishes,  is  shown  for 
reference  . 

The  general  observation  to  be  made  from  Figures  2  through  5  is  that 
the  approximate  dispersion  relation  provides  a  highly  satisfactory 
representation  to  the  exact  dispersion  relation  for  the  submerged  jet 
throughout  wave  number  space.  The  representation  is  best  for  the  weakest 
and  widest  jets  and  degrades  with  increasing  jet  strength  and  decreasing 
jet  width  . 

A  second  observation  to  be  made  from  these  figures  is  that  a  totally 
submerged  current  field  can  appreciably  modify  the  dispersion  relation 
for  a  surface  wave  vis-a-vis  its  dispersion  relation  in  a  region  where  the 
current  field  vanishes  .  This  fact  could  have  important  implications  for 
the  understanding  and  interpretation  of  oceanic  features  sensed  by  radar . 

CONCLUSIONS 

An  approach  has  been  developed  for  the  analysis  of  the  modifications 
of  surface  waves  by  depth-varying  current  fields  .  The  approach  is  based 
on  an  approximate  dispersion  relation  for  wave-current  interactions  derived 
from  the  governing  equations  of  the  problem  (the  inviscid  Orr-Sommerfeld 
equations)  coupled  with  the  wave  kineraat tr/wave  action  formulation  of 
surface  wave  propagation.  The  wave  kinematic /wave  action  formulation 
arises  purely  from  kinematic  and  dynamic  considerations;  hence,  the 
accuracy  of  the  approach  depends  entirely  on  the  accuracy  of  the  approx¬ 
imate  dispersion  relation  with  respect  to  the  exact  dispersion  rela¬ 
tion.  We  have  demonstrated  here  and  elsewhere  (Skop  (5))  that  the 
approximate  dispersion  relation  provides  a  highly  satisfactory  repre- 


sentation  to  the  exact  dispersion  relation  for  a  variety  of  depth- 


varvtng  current  fields  ranging  from  simple  to  complex  In  structure. 

The  approach  Is  particularly  useful  In  that  It  focuses  directly 
on  the  wave  parameters  of  Interest  (the  amplitude,  frequency,  direction, 
and  wavelength  of  the  wave)  and  eliminates  the  requirement  to  solve  the 
Invlsctd  Orr-Sommerfeld  equations  to  derive  these  parameters. 
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APPENDIX  —  DISPERSION  RELATION  FOR  THE  SUBMERGED  JET 


With  U2  3  0  and  =  k,  the  effective  velocity  U  defined  by  equation 
(11)  becomes  U  ■  •  Noting  from  equation  (26a)  that  3^U^/3z^  -  0, 

equation  (12a)  for  the  vertical  velocity  component  w  becomes 


92w  _ 

-  k2w  =  0 

3z2 


(Al) 


Since  U3  =  0  and  3U^/3z  =  0  at  the  free  surface,  the  boundary  condition 
there  is,  from  equation  (12b)  with  y  =  0, 

~  3w 

—  =  gw  at  z  =*  0  (A2) 

3z 


Discontinuities  in  the  velocity  gradient  of  the  submerged  jet  occur  at 
z  *  -  ad,  -d,  and  -(2  -  a)d  .  Across  these  discontinuities,  the 
vertical  velocity  and  the  pressure  must  be  continuous.  Hence,  we  have 

3w  3UX 

w  and  (U^  -  c) - w  (A3) 

3z  3z 

must  be  continuous  at  z  *  -  ad,  -  d,  and  -(2  -  a)d . 

We  seek  solutions  for  equation  (Al)  as 

„  „kz  .  ■  „  -k  z 
ale  +  “le 

a2ek(z+od)  +  ^^(z+ad) 

a3ek(z+d)  +  b3e“k(z+d) 
a4ektz+(2-a)d, 


0  >  z  >  -ad 
-ad  >  z  >  -ad 
-d  >  z  >  -(2  -  a)d 
-(2  -  a)d  >  z 


(A4) 


where  a 3  through  a^  and  b^  through  bj  are  constants  of  integration.  Apply- 
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lag  the  boundary  and  continuity  conditions  specified  by  equations  (A2) 
and  (A3),  recalling  that  c  =  u>/k ,  and  introducing  the  dimensionless 
parameters  defined  by  equation  (28),  we  find 


0^(a^  -  b^)  =  k(s^  "t  b^) 

(A  5a) 

a^e  a<  +  b^eaK:  =  a2  +  b2 

(A  5b) 

F 

0(aie-a<  -  b1eaK)  =  JJ(a2  -  b2)  + - (a2  +  b2) 

(1-a) 

(A5c) 

a2e~(  1-<X) K  +  b2e(1-a>K  =  a3  +  b3 

(A5d) 

(Fk  +  0)[a2e-(1-a)K  -  b2e(1“a)Kl 

+  ----  la2e“(1'a),c  +  b2e(1"a)K] 
a-a) 

F 

=  (Fk  +  tl)(a3  -  bj)  - - (a3  +  b3) 

(1-ot) 

(A5e) 

a3e-(1-a>*  +  b3e<1-a),c  =  a4 

(A5f  ) 

n[a,e“(1“a),e  -  b3e(1"a>K]  - 

(1  -  a) 

+  b3  e(1-°t>K]  =  0  a4 

(A5g) 

Elimination  of  the  constants  of  integration  from  equations 

(A5)  leads  to 

the  exact  dispersion  relation  given  by  equation  (29). 


20 


